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ABSTRACT 

We use high resolution (~ 0.1") F814W ACS images from the Hubble Space Telescope ACS Treasury 
survey of the Coma cluster at z ~ 0.02 to study bars in massive disk galaxies (SOs), as well as low-mass 
dwarf galaxies in the core of the Coma cluster, the densest environment in the nearby Universe. Our study 
helps to constrain the evolution of bars and disks in dense environments and provides a comparison point 
for studies in lower density environments and at higher redshifts. Our results are: (1) We characterize 
the fraction and properties of bars in a sample of 32 bright (My < —18, M* > 1O 9 - 5 M ) SO galaxies, 
which dominate the population of massive disk galaxies in the Coma core. We find that the measurement 
of a bar fraction among SO galaxies must be handled with special care due to the difficulty in separating 
unbarred SOs from ellipticals, and the potential dilution of the bar signature by light from a relatively 
large, bright bulge. The results depend sensitively on the method used: the bar fraction for bright SOs 
in the Coma core is 50± 11%, 65± 11%, and 60± 11% based on three methods of bar detection, namely 
strict ellipse fit criteria, relaxed ellipse fit criteria, and visual classification. (2) We compare the SO bar 
fraction across different environments (the Coma core, Abell 901/902, and Virgo) adopting the critical 
step of using matched samples and matched methods in order to ensure robust comparisons. We find that 
the bar fraction among bright SO galaxies does not show a statistically significant variation (within the 
error bars of ±11%) across environments which span two orders of magnitude in galaxy number density 
(n ~ 300-10,000 gal/Mpc 3 ), and include rich and poor clusters, such as the core of Coma, the Abell 
901/902 cluster, and Virgo. We speculate that the bar fraction among SOs is not significantly enhanced 
in rich clusters compared to low density environments due to two reasons. Firstly, SOs in rich clusters 
are less prone to bar instabilities as they are dynamically heated by harassment and are gas poor as a 
result of ram pressure stripping and accelerated star formation. Secondly, high-speed encounters in rich 
clusters may be less effective than slow, strong encounters in inducing bars. (3) We also take advantage 
of the high resolution of the ACS (~ 50 pc) to analyze a sample of 333 faint (My > —18) dwarf galaxies 
in the Coma core. Using visual inspection of unsharp-masked images, we find only 13 galaxies with bar 
and/or spiral structure. An additional eight galaxies show evidence for an inclined disk. The paucity of 
disk structures in Coma dwarfs suggests that either disks are not common in these galaxies, or that any 
disks present are too hot to develop instabilities. 
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rals (Weinzirl et al. 2009; Kormendy & Fisher 2008; Lau- 
rikainen et al. 2007) since z < 1. Other processes, such as 
minor mergers (Jogee et al. 2009; Weinzirl et al. 2009), 
or internal secular processes (Laurikainen et al. 2007; Ko- 
rmendy & Fisher 2008) are increasingly invoked. 

The most efficient internal driver of evolution in disk 
galaxies are stellar bars. They effectively redistribute an- 
gular momentum in the disk and dark matter (DM) halo, 
and drive large gas inflows to the central regions of galax- 
ies. The resulting central dense gas concentrations ignite 
powerful starbursts (Schwarz 1981; Sakamoto 1999; Jogee 
1999; Jogee et al. 2005; Sheth et al. 2005). In this way, 
bars are thought to build disky central components known 
as pseudobulges (Hohl 1975; Kormendy 1979; Combes & 
Sanders 1981; Combes 1990; Kormendy 1993; Jogee 1999, 
Jogee et al. 2005). 

We describe below some of the recent progress made in 
exploring the bar fraction (defined as the fraction of disk 
galaxies hosting a large-scale bar) in terms of methodol- 
ogy, dependence on Hubble type, redshift, and environ- 
ment. For a long time, statistics on the optical bar frac- 
tion in the local universe came from visual classification 
of the galaxies in the Third Reference Catalog of Bright 
Galaxies (RC3; deVaucouleurs et al. 1991). In the RC3, 
over all disk galaxies (SO-Im), the visual optical bar frac- 
tion is ~ 30% for strong bars (SB), and ~ 30% for weak 
bars (SAB), giving ~ 60% overall. The RC3 optical bar 
fraction suffers from two limitations: firstly, it denotes the 
optical bar fraction averaged over a broad range of Hub- 
ble types, and secondly, it is based on visual classification 
to identify and characterize bars, thus giving no quantita- 
tive measurements of their properties (such as size, shape, 
strength). In fact, the former limitation has persisted in 
many recent studies which focus on the average bar frac- 
tion over disks of all Hubble types (Eskridge et al. 2000; 
Knapen et al. 2000; Hunt & Malkan 1999; Mulchaey & 
Regan 1997; Laine et al. 2002; Laurikainen et al. 2004; 
Block et al. 2004; although see Odewahn 1996). 

Recent studies have made important headway on both 
fronts. Firstly, quantitative methods such as ellipse fit- 
ting (e.g., Wozniak et al. 1995, Friedli et al. 1996; Re- 
gan & Elmegreen 1997; Mulchaey & Regan 1997; Jogee et 
al. 1999, 2002, 2004; Knapen et al. 2000; Laine et al. 2002; 
Sheth et al. 2003, 2008; Elmegreen et al. 2004; Menendez- 
Delmestre et al. 2007; Marinova & Jogee et al. 2007, 
hereafter MJ07; Aguerri et al. 2009), and bulge-bar-disk 
decomposition (e.g., Laurikainen et al. 2005; Reese et 
al. 2007, Weinzirl et al. 2009, Gadotti 2011) are being used 
to reduce the level of subjectivity. Secondly, evidence is 
mounting that the bar fraction varies across galaxies of 
different Hubble types (Barazza, Jogee, & Marinova 2008, 
hereafter BJM08; Aguerri et al. 2009; Marinova et al. 2009, 
hereafter M09) and depends non-monotonically on the 
host galaxy properties, such as bulge-to-total ratio (B/T), 
luminosity, stellar mass, and color, at a range of redshifts 
and environments (M09; Barazza et al. 2009; Weinzirl et 
al. 2009; Laurikainen et al. 2009; Nair & Abraham 2010a; 
Cameron et al. 2010; Gadotti 2011). 

Most of the results described above have focused on field 
galaxies at low redshifts, while much less is known about 
barred disks at higher redshifts (Abraham et al. 1999; 
Sheth et al. 2003; Jogee et al. 2004; Sheth et al. 2008; 



Cameron et al. 2010), and in dense environments. Dense 
environments such as clusters can be an important lab- 
oratory for studying the co-evolution of bars and their 
host disks, as there are many physical processes, which are 
unique to such environments and impact disks and bars. 
For example, we can study how the cumulative effect of 
frequent weak galaxy encounters ('harassment'; Moore et 
al. 1996), the effect of the intra-cluster medium (ICM) on 
galaxies (Gunn & Gott 1972; Larson et al. 1980; Quilis 
et al. 2000) , and the influence of the cluster potential as a 
whole (e.g., Byrd & Valtonen 1990; Gnedin 2003) affect the 
properties of bars and disks. Theoretical studies exploring 
these processes give conflicting predictions. For example, 
although the numerous tidal encounters in a dense cluster 
core can induce bars in unbarred disks in the case of a 
prograde encounter, they can also have little or no effect 
in terms of inducing a bar (or affecting the strength of an 
already existing bar) in the case of retrograde encounters 
(Gerin, Combes, & Athanassoula 1990; Romano-Diaz et 
al. 2008; Aguerri & Gonzalez- Garcia 2009). Determining 
the fraction and properties of barred galaxies in such ex- 
treme environments and comparing them to those of galax- 
ies in environments of varying densities, can give clues 
to the outstanding problems in understanding the forma- 
tion and growth of bar instabilities in disks, and therefore 
to understanding disk evolution. In addition, identifying 
barred disks in cluster environments can provide a lower 
limit to the fraction of disk systems in clusters (e.g., M09; 
Mendez-Abreu et al. 2010). 

Quantitative results addressing this issue are only start- 
ing to emerge. Several recent studies of non-dwarf galax- 
ies at a range of redshifts from 0.01 < z < 0.4 (Barazza 
et al. 2009; Aguerri et al. 2009; M09), find that the opti- 
cal bar fraction shows at most a modest variation (±10%) 
between field and intermediate density environments (e.g., 
moderately rich clusters) . But what happens at really high 
densities (i.e., dense cluster cores)? Some studies (Thomp- 
son 1981; Andersen 1996; Barazza et al. 2009) have sug- 
gested that, within a galaxy cluster, the bar fraction may 
be higher in the dense core regions than the outskirts. 
However, this has remained an open question for cluster 
cores due to issues such as limited number statistics, pro- 
jection effects, poor quality or inhomogeneous data, and 
uncertainties in cluster membership. We make progress in 
two ways with this study. Firstly, we are able to establish 
cluster membership for all galaxies in our bright SO sample 
(see § [3]) using spectroscopic redshifts. Although we are 
still limited by line-of-sight projection effects from within 
the Coma cluster itself, the number density of galaxies in 
the core of Coma is a factor of 10-100 times higher than 
at larger cluster radii, so we expect the impact of such 
linc-of-sight contamination to be low. Secondly, the high- 
resolution HST ACS data is vital in identifying bars in 
SOs, especially in the types of cases we discuss in Sec- 
tion 13.2.21 namely, where a relatively large bulge, in com- 
bination with a relatively short bar and projection effects 
combine to make bar detection difficult in low-resolution 
ground-based images. 

In this paper, we provide baseline results for the densest 
environment in the local Universe: the central regions of 
Coma. We conduct two explorations using two different 
samples. In the first part of this paper (§ [3]), we focus on 
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the sample of bright (M\ 



< 



-18) non-dwarf SO galaxies, 



which comprise 94% of our sample of bright disk galax- 
ies in the Coma core. We characterize their disk and bar 
properties and present results on the optical bar fraction 
for bright SO galaxies from ellipse fitting and visual classifi- 
cation. We compare with the results of other studies of SO 
galaxies in Coma, and less dense clusters (Abell 901/902 
and Virgo). 

In the second part of this paper (§ 2]), we take advan- 
tage of the exquisite 50 pc resolution of the ACS images 
to search for bars and other disk features (spiral arms, in- 
clined disks) in the faint, dwarf (My > —18) galaxies of 
the Coma cluster core. The prevalence or paucity of bar 
structures in early-type dwarf galaxies can not only pro- 
vide clues on the evolutionary history of these dwarf sys- 
tems, but also has implications for the conditions necessary 
for bar formation and growth in galaxies. Are bar/spiral 
arm instabilities commonplace for dwarf galaxies? Some 
early-type dwarfs in clusters are believed to have origi- 
nated from late-type spirals or dwarf irregulars (e.g., Lin 
& Faber 1983; Kormendy 1985) that have been stripped of 
their gas by external processes. This is especially true in 
a cluster setting where environmental processes like ram- 
pressure stripping and harassment are commonplace. If 
the progenitor galaxy hosted a stellar bar or spiral arms, 
these stellar features can persist even after the galaxy's gas 
has been processed by the cluster environment. In fact, a 
number of early-type dwarf galaxies in Virgo have been 
observed to host stellar disk features such as a lens, a bar, 
or spiral arms (e.g., Sandage & Binggeli 1984; Binggeli & 
Cameron 1991; Jerjen et al. 2000; Barazza et al. 2002; 
Lisker et al. 2006; Lisker et al. 2007; Lisker & Fuchs 2009) 
thus supporting the above scenario. Graham et al. (2003) 
discovered the first two early-type dwarf galaxies with 
spiral-like structure in the Coma cluster using unsharp 
masking, concluding that such galaxies may provide the 
'missing link' in the evolution of faint spiral galaxies into 
dwarf spheroidals due to cluster processes. 

The layout of this paper is as follows: In § we de- 
scribe our dataset and sample selection. Section [3] deals 
with our bright (My < —18) Coma core sample. In this 
section we outline our methods for identifying bright SO 
galaxies (§ 13. 1[) . our methods for identifying bars in these 
galaxies (§ 13.21 and 13. 3[) . our results for SO galaxies in the 
Coma core (§ 13. 4K l3.7p . and the implications of these re- 
sults for bar and disk formation and evolution in dense 
environments (§ 13. 8p . In §[51 we present our investigation 
of the faint, dwarf galaxies in our Coma core sample. We 
describe our methods for finding disk features (bar, spiral 
arms, inclined disk) in these galaxies (§ 14. 2p , as well as 
our results and discussion for dwarf galaxies in § 14.31 We 
summarize all of our results in § [5j 

2. DATA AND SELECTION OF A CLUSTER SAMPLE 



Our data come from the Hubble Space Telescope (HST) 
Advanced Camera for Surveys (ACS) Treasury survey of 
the Coma cluster at z ~ 0.023 (Carter et al. 2008). Orig- 
inally designed to cover a large area of the core and infall 
region of Coma, the survey remains only ~ 28% complete 
because of the failure of ACS in 2007. Nevertheless, the 

16 We assume in this paper a flat cosmology with Qm = 1 — Qjy 

17 The transformation equation tables can be found at |http:/ /w 



available data span 274 arcmin 2 , where approximately 75% 
of the data are within 0.5 Mpc of the cluster center. The 
data cover approximately 70% of the core region of Coma 
(assuming R COIC ~ 0.2 Mpc) and are therefore representa- 
tive of the core and the immediate surroundings, namely 
the region where the galaxy number density is ~ 10,000 
galaxies Mpc~ 3 before it drops sharply as a function of dis- 
tance from the cluster center (The & White 1986). This 
dataset contains thousands of sources down to a limiting 
magnitude of I = 26.8 mag in F814W (AB mag). The ACS 
point-spread function (PSF) is ~ 0.1", which corresponds 
to ~ 50 pc at the distance of the Coma cluster, assum- 
ing a distance of D = 100 MpcQ SExtractor source 
catalogs are available as a part of the Coma survey data 
releases. The second data release (DR2) is described in 
detail in Hammer et al. (2010). Throughout the paper, 
Mi(g 14 ) as well as SDSS g and r magnitudes are given in 
the AB system, while B and V magnitudes are in the Vega 
system. 

We first use the eyeball catalog of Trentham et al. 
(in preparation) to select cluster members. In this cat- 
alog, galaxies are visually assigned a cluster membership 
class from to 4. Galaxies with membership class are 
spectroscopically confirmed members, while galaxies with 
class 4 are visually deemed to be likely background ob- 
jects. The intermediate membership classes from 1 (very 
probable cluster member) to 3 (plausible cluster member), 
are assigned based on a visual estimation taking into ac- 
count both surface brightness and morphology (Trentham 
et al. in prep). We select objects with apparent mag- 
nitude mi(8i4) < 24 (AB mag) and membership class 
to 3, resulting in a sample of 469 cluster galaxies. For 
these galaxies, 41% are spectroscopically confirmed mem- 
bers (class 0), while 7%, 27%, and 25% have membership 
classes 1 to 3, respectively. 

We derive B and V magnitudes (in Vega mag) for the 
Coma galaxies using SDSS g and r (in AB mag). For the 
bright sample, we use the B, V, g, and r magnitudes (in- 
stead of the ACS F814W) for ease of comparison to other 
studies. In addition, Hammer et al. (2010) find that for 
bright galaxies (M I ( 814 ) < 17 AB mag) in the Coma sur- 
vey, it is more reliable to use the SDSS rather than ACS 
magnitudes, as the latter may be unreliable for some galax- 
ies with large, diffuse stellar halos. We use the following 
transformations from Jester et al. (2005) to convert the 
SDSS g and r (AB) to B and V (Vega^J: 

5 = a + 0.39x fo- 70+0.21 



g + U.ciy x (g - 
V = g -0.59 x(g-r) -0.01. 



(1) 
(2) 

a distance 



We calculate absolute magnitudes assuming 
modulus of 35.0 (Carter et al. 2008). 

In this paper, we explore the optical bar fraction in two 
regimes: bright, non-dwarf SO disks (§[3]) and faint (dwarf) 
galaxies (§[3|. To separate these two regimes, we apply a 
magnitude cut of Mj( 81 4) < —18.5, roughly equivalent to 
My < — 18 or Mb < — 17 for our sample. We choose a cut 
at My ~ — 18 because it tends to separate well the regimes 
where normal and dwarf galaxies dominate on the lumi- 
nosity functions of clusters (Binggeli et al. 1988; Trentham 
1988; Trentham & Hodgkin 2002; Mobasher et al. 2003). 

0.3 and H =70 km s" 1 Mpc" 1 . 

'.sdss.org/dr7/algorithms/sdssUBVRITransform.html 
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Luminosity cuts at magnitudes My ~ — 18 are often used 
to separate dwarf and non-dwarf galaxies in the literature 
(Matkovic & Guzman 2005; Aguerri et al. 2005; van Zee 
2001; Barazza et al. 2006). This cut gives 52 galaxies 
brighter than, and 417 galaxies fainter than My = —18. 
From the 52 bright galaxies, we discard four galaxies (two 
are a close/merging pair and two galaxies are partially off 
the edge of a tile) bringing our initial bright sample to 48 
galaxies. We discuss the methods for selecting SO galaxies 
from this bright sample in § 13.11 

In Fig. [IJt, we show the absolute magnitude Mt( 81 4) dis- 
tribution for the non-dwarf (bright) and dwarf Coma core 
samples. Fig. [TJ) shows the peak surface brightness /xo 
(from the SExtractor source catalogs in DR2) vs. ab- 
solute magnitude Mi(gi4) distribution of the bright, non- 
dwarf and dwarf galaxies. 

3. BARS IN BRIGHT SO GALAXIES IN THE CENTRAL 
REGION OF THE COMA CLUSTER 

3.1. Identifying Bright SO Galaxies 

Due to the fact that bars are inherently a disk phe- 
nomenon, the bar fraction is traditionally defined as 

, _ iV ba 

/bar 



iVbar + N unh 

ar 

where A bar and -/V un b ar represent the number of barred 
and unbarred disk galaxies, respectively. Therefore, from 
the bright sample of 48 galaxies, we need a sample of disk 
galaxies (e.g., SO-Im) for analysis. As discussed in § [TJ 
recent work has shown that a bar fraction averaged over a 
wide range in Hubble types gives only limited information. 
The bar fraction is a strong function of galaxy properties, 
such as B/T, luminosity, stellar mass, and color. Because 
our sample is too small to split into fine bins by morpholog- 
ical type, and because most (94%) of bright disk galaxies in 
our Coma core sample are SOs (see below), our analysis of 
bright galaxies in this paper focuses on SOs only. The goal 
of our study is to provide the bar fraction for the densest 
low redshift (z ~ 0.02) environment and to serve as a com- 
parison point for studies of barred SO galaxies in field and 
intermediate-density environments at different redshifts. 

Starting with the bright sample of 48 galaxies, we use 
visual classification to separate the galaxies into ellipti- 
cals, SOs, visually ambiguous E/S0, and spirals. We note 
that our visually-identified class of 'SO' galaxies includes 
all Hubble type SO sub-types from S0~ to SO/a (numer- 
ical T-types -3 to 0). It is fairly easy to separate SOs 
and ellipticals visually when the SOs host bars. However, 
unbarred SOs are harder to separate from ellipticals since 
an unbarred SO hosts a disk, which is effectively feature- 
less and devoid of tell-tale disk signatures, such as a bar 
or spiral arms. We find a group of 10 bright galaxies 
that are visually ambiguous E/SOs. For each of these 10 
galaxies, we perform multiple component structural de- 
composition with the GALFIT code (Peng et al. 2002) by 
fitting the two-dimensional (2D) light distribution taking 
into account the PSF, following the procedure described 
in Weinzirl et al. (2009). In brief, we fitted each galaxy 
with three models: a single-component Sersic model, a 
two-component ('bulge'+disk) model, and a three compo- 
nent ('bulge'+disk+bar) model. The disk is represented 



by an exponential (Sersic n = 1) model, while the 'bulge' 
is a Sersic model with a free-floating Sersic index n. The 
bar is typically associated with a Sersic model of low n. 
If needed, a point source component was added to repre- 
sent the point sources, which are common in the center of 
these galaxies. For each model, GALFIT finds the opti- 
mum solution using the Levenberg-Marquardt algorithm. 
The goodness of fit is determined iteratively by calculat- 
ing x 2 - GALFIT continues to adjust the model parame- 
ters until the gradient 5x 2 /x 2 is ver Y small (e.g., 10~ 4 ) for 
10 continuous iterations. Out of the three models (Sersic, 
'bulge'+disk, 'bulge'+disk+bar), the best one was selected 
by considering a number of factors, including the \ 2 va l~ 
ues, the strength and spatial distribution of the residuals, 
and the physical viability of output parameters (e.g., effec- 
tive radii of the Sersic and bulge component; scale length, 
axial ratio, and position angle of the disk). These factors 
were used to decide whether the galaxy is likely an ellipti- 
cal or an SO. Examples of the data, model, and residuals 
for three representative galaxies classified as E, SO/a, and 
E/S0 after decomposition are shown in Figured 

Out of the 10 visually ambiguous E/S0 galaxies, we find 
from multi-component decompositions that one is an el- 
liptical, eight are SO or SO/a, and one is still ambiguous 
E/S0. The ambiguous galaxy may be a disk galaxy with 
complex structure or an elliptical with an inner debris disk. 
In § 13.41 we use this galaxy to estimate the uncertainty in 
the optical bar fraction by calculating /bar (see Equation 
1) for the two cases where this galaxy is either included or 
excluded in the number of unbarred disks (iV un b a r)- 

The absolute magnitude AiT( 814 ) distribution of the 
bright Coma core sample is shown in Fig. [3^. The final 
morphological breakdown of our bright sample (13 ellip- 
ticals, 1 ambiguous E/S0, 32 SOs, and 2 spirals) is shown 
in Fig[3}D. It is clear that SOs dominate among the bright 
disk galaxies in our Coma core sample, which is expected 
for the central regions of a dense cluster. We find a ratio of 
E : SO : Sp of 28% : 68% : 4%. This is at the extreme end 
of the morphology-density relation found in dense environ- 
ments by Dressier (1980). Fig. [3J; shows the distribution 
of stellar mass for the SO disk sample as well as all bright 



(M 



< 



-18) galaxies. Stellar masses are calculated using 



the relations from Borch et al. (2006)0 assuming a Kroupa 
et al. (1993) initial mass function: 



M 



x 10 



-0. 628+1. 305(B-V)-0.1 



where 



''lu 



10 



-0.4(V-4.82) 



(4) 



(5) 



Galaxies in our SO disk sample have stellar masses between 
lO 9 ' 5 and 10 n M Q . Fig. Hi shows a g - r color vs. M r 
magnitude diagram. Almost all ellipticals and most disk 
galaxies fall on the red sequence. We overplot the relation 
from Blanton et al. (2005a) for the break between the red 
sequence and blue cloud using the equation (modified with 
an offset of -0.77 for h = 0.7) 

( g - r )= 0.65 - 0.03(M r - 0.77 + 20). (6) 

We also plot a subsample of the dwarf galaxies (My > 
-18) for which SDSS data are available (~ 30%). 



18 The Kroupa IMF offset term is reported as -0.15 in Bell et al. (2003). However this value was calculated assuming unrealistic conditions 
(Bell, E., private communication). The correct value of -0.1 was recalculated and reported in Borch et al. (2006). 
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We show examples from our final bright SO sample of 32 
galaxies in Figure |4j We note that all 32 SOs in the bright 
sample are spectroscopically confirmed cluster members. 

3.2. Identification of bars in SOs via ellipse fits 

Ellipse fitting is our primary method of detecting bars 
in the bright SO sample (e.g., Wozniak et al. 1995, Fricdli 
et al. 1996; Regan et al. 1997; Mulchaey & Regan 1997; 
Jogee et al. 1999, 2002, 2004; Knapen et al. 2000; Laine 
et al. 2002; Sheth et al. 2003, 2008; Elmegreen et al. 2004; 
Menendez-Delmestre et al. 2007; MJ07; Aguerri et 
al. 2009). To detect bars through ellipse fitting we use 
the standard IRAF task ELLIPSE in conjunction with 
an adaptive wrapper (Jogee et al. 2004), which runs 
ELLIPSE iteratively on each galaxy until the best fit 
is found or up to a maximum number of times specified 
by the user. Ellipses are fit to the galaxy isophotes out 
to a maximum distance (adisk) where the brightness of the 
isophotes reaches the noise level. We note that the value of 
Odisk depends on the depth of the image, as adisk will reach 
larger values for deeper images. However for the purpose 
of bar detection, it is only necessary for the radial profile 
to extend beyond the bar into the more circular region of 
the disk. We typically set the maximum allowed iterations 
to 300, however for most galaxies a good fit is achieved in 
only a few iterations. A good fit is one where an ellipse can 
be fitted at every isophote out to adisk- Residuals charac- 
terizing how well each isophote is fitted by its correspond- 
ing ellipse are given by the harmonic amplitudes A3, B3, 
A4, and B4 (e.g., Carter 1978; Jedrzejewski 1987; Carter 
1987). For our barred galaxies, we find typical amplitudes 
of 5-10%. For a detailed discussion on the advantages and 
drawbacks of using ellipse fitting to characterize bars we 
refer the reader to MJ07. 

Once the galaxies are fitted, we use an interactive visual- 
ization tool to display the overlays of the fit on the galaxy 
image, as well as the radial profiles of surface brightness, 
cllipticity (e), and position angle (PA). Using the radial 
profiles of ellipticity (e) and PA, we classify the galaxies 
as 'highly inclined', 'barred', or 'unbarred'. We discuss 
these classes in more detail below. 

3.2.1. Detecting and removing highly inclined galaxies 

In studies of bars, it is conventional to exclude highly 
inclined galaxies as the large inclination precludes accu- 
rate structural classification, making it particularly diffi- 
cult to identify systems as barred or unbarred. We use two 
ways to identify highly inclined galaxies in the bright SO 
sample. The first is via the ellipse fit criteria, where the 
observed outermost disk isophote (at adisk) has e > 0.5, 
corresponding to i > 60°. We find eight galaxies that 
fit this criterion. This method works well for spirals of 
intermediate to late Hubble types, but does not capture 
all highly inclined disks for SOs because for some edge-on 
or highly inclined SOs, a rounder, thickened outer stellar 
component can sometimes dilute the outermost isophote 
so that the outermost ellipticity is below 0.5 although the 
galaxy is highly inclined. 

Therefore, our second method is to visually identify 
highly-inclined SOs. We visually identify these systems 
using the criteria that a thinner, high-surface brightness, 
highly-inclined disk appears embedded in a thick, diffuse 



stellar component, which could be a mix of thick disk, 
bulge, and bar stars. A typical example of one of these 
galaxies is shown in Fig. \5\ We encounter four such cases, 
where the galaxy appears visually to be close to edge-on, 
and is classified as highly inclined. We distinguish these 
cases from a face-on galaxy with a bar because (1) the thin, 
highly-inclined disk and the thicker stellar component are 
always oriented along the same position angle, (2) the thick 
outer stellar component in these highly-inclined SO galax- 
ies appears much fainter and more diffuse than a face-on 
disk, and (3) in three out of these four galaxies, a box or 
X-shaped bulge is present, suggesting that the galaxy is 
seen edge-on (e.g., Athanassoula 2005). We therefore ex- 
clude from further analysis the 12 highly- inclined systems 
that we find in the sample. 

3.2.2. Detecting barred galaxies 

Traditionally, when ellipse fits are used to identify bars, 
a galaxy is classified as barred if the radial profiles of the el- 
lipticity and PA fulfill the following requirements: (1) the 
e rises to a global maximum, ebar > 0.25, while the PA 
remains relatively constant (within ±10°), and (2) the e 
drops by at least 0.1 and the PA changes by more than 
10° at the transition between the bar and disk region. An 
example of a barred SO galaxy in our sample that meets 
the traditional criteria is shown in Fig. [BJ 

These criteria will identify primary stellar bars in the 
vast majority of spirals, particularly those of intermedi- 
ate to late Hubble types (Sb-Sm). However, they can 
marginally fail in some galaxies due to a rare set of circum- 
stances, which we describe in detail below. These circum- 
stances are particularly likely to occur in SOs with large 
bulgc-to-disk ratios. 

In some SOs, a combination of structural parameters and 
viewing angle causes the observed (i.e., not deprojected) 
maximum bar ellipticity (ebar! typically measured from 
the isophote crossing the end of the bar) to become a local 
maximum of the ellipticity radial profile rather than the 
global maximum. In these rare cases, the ellipticity (edisk) 
of the outer disk becomes the global maximum in the ra- 
dial profile of ellipticity. This can happen in the case of 
a barred galaxy, where all or most of the following condi- 
tions are satisfied: (i) The galaxy has a moderate to large 
inclination (e.g., i > 50°). This causes the outer circu- 
lar disk of the galaxy to appear elongated along the line 
of nodes (LON) in the projected image of the galaxy on 
the sky, leading to a higher measured edisk- (h) A large 
fraction of the length of the bar lies within a fairly axisym- 
metric bulge, which is much more luminous than the bar. 
In this case, the bulge light dilutes the ellipticity of the bar 
by 'circularizing' the isophote crossing the bar end, thus 
causing ebar, measured from this isophote to be signifi- 
cantly lower than the true ellipticity of the bar. (hi) The 
bar major axis has a large offset (A6*) with respect to the 
LON, such that projection effects make the disk appear 
more elongated, while the bar appears more round. The 
most extreme example occurs when the bar is perpendic- 
ular to the LON (i.e., A8 = 90°). Such situations can 
potentially cause the observed ellipticity of the disk to ex- 
ceed that of the bar. Thus, a combination of factors (i) to 
(iii) can cause the measured ebar to fall below edisk so that 
ebar is a local maximum in the e radial profiles. In this 
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case, the bar can still be identified through ellipse-fits if 
the traditional criterion that the measured maximum bar 
cllipticity ebar must be a global maximum is relaxed, and 
a local maximum be deemed acceptable. 

In the case of barred SOs, the conditions (i) to (hi) can be 
satisfied in a larger fraction of galaxies than for a sample of 
barred intermediate-to-late Hubble type (Sb-Sm) spirals 
due to the following reasons. Many barred cluster SOs host 
bulges that are bright, have large bulge-to-disk light ratios, 
and encompass a large fraction of the length of the bar. In- 
deed, among our sample of 20 moderately-inclined cluster 
SOs, we find three such cases and an example is shown in 
Fig. [7] For this reason, we quote two bar fractions derived 
through ellipse fits: the first bar fraction (/bar.Es), where 
we use the strict criteria (1) and (2) above, and the second 
bar fraction (/bar,ER) where for galaxies satisfying (i) to 
(hi), we relax the criterion that the maximum bar ellip- 
ticity must be a global maximum (however we still require 
it to rise above e = 0.25). We note that all bars iden- 
tified with the strict ellipse-fitting criteria ('ES') are also 
picked up under the relaxed ellipse-fitting criteria ('ER'). 
We further note that, for the galaxies where the bar is 
detected only through the relaxed criteria, if the radial 
profiles of the ellipticity and PA are deprojected (§ 13. 7p . 
the bar ellipticity ebar does become the global maximum. 
However, because many large studies of bars do not de- 
project the radial profiles (e.g., M09), we use the observed 
radial profiles to detect the bars as described above for 
ease of comparison. 

3.3. Identification of bars in SOs via visual classification 

In addition to ellipse fitting, we also present the optical 
bar fraction for bright SOs in Coma from visual classifi- 
cation performed by I.M., S.J., and P.E. This facilitates 
comparison to other work where bars are identified visu- 
ally (§& 

A galaxy is classified as 'barred' through visual classi- 
fication if it has a significant elongated feature extending 
from the center of the disk with an axial ratio (estimated 
from the image with DS9) < 0.7 and a PA that differs 
from the PA of the outer disk by at least 10°. A galaxy 
is classified as 'unbarred' if there is no elongated structure 
present that fits the above criteria. 

All of the bright barred SOs we identify in the Coma 
sample through ellipse fitting and visual classification are 
listed in Table [TJ The methods through which the bar is 
detected are shown in column (5). 

3.4. Optical SO bar fraction in the central region of the 

Coma cluster 

The optical bar fractions for our sample of bright SO 
galaxies in the Coma cluster core are presented in Tabled] 

Using the strict ellipse fitting criteria (§ I3.2.2[) . we 
find that the optical bar fraction for the bright SOs is 
/bar : ES = 50 ± 11% (10/20). Using the relaxed ellipse 
fitting criteria, we find /bar,ER = 65 ± 11% (13/20). Vi- 
sual classification gives an optical bar fraction of 60 ± 11% 
(12/20). All errors are binomial errors. The barred SO 
galaxies identified through ellipse fitting and visual classi- 
fication are shown in Fig. [5] 

19 This correction is used by T81 to account for the effects of foregroum 
of spectroscopic data. Since all bright SOs in our Coma core sample 
corrected bar fraction from T81. 



To correctly derive the bar fraction for SOs in clusters, 
we need to accurately estimate the number of unbarred SOs 
(-^unbar in Eq. 1). As SOs are devoid of typical disk fea- 
tures such as spiral arms, star- forming rings, etc., it is par- 
ticularly challenging to visually identify all unbarred SOs 
and separate them from ellipticals. Therefore, in § 13.11 we 
identified SOs through both visual classification and two- 
dimensional structural decomposition of the images into 
single-component Sersic models, bulge+disk models, and 
bulge+disk+bar models. We found 13 Es, 32 SOs, and one 
E/S0 case, which still remains ambiguous even after de- 
composition. This ambiguous E/S0 case is not included in 
the optical bar fraction in Table [5] If it is included as an 
unbarred SO in our analysis, the optical bar fractions fall 
to: /bar.ES = 48± 11%, / bar , E R = 62± 11%, and 57± 11% 
from visual classification. We therefore estimate that the 
uncertainties associated with determining the number of 
unbarred SOs can lead us to overestimate the optical bar 
fraction by only a small factor of ~ 1.05. 

We note that in the field, the bar fraction is lower in 
the optical than in the NIR by factor of 1.3 (Eskridge et 
al. 2000; M07) for intermediate (Sab-Sc) Hubble types due 
to obscuration by gas, dust, and star formation. However 
in SOs, where there is little gas and dust on large scales, we 
don't expect the difference between the optical and NIR 
bar fractions to be significant. 

3.5. SO bar fraction across different environments 

Due to the fact that different bar detection methods can 
yield different bar fraction results, it is important to com- 
pare studies using the same methods for consistency. In 
addition, as discussed in § [TJ the bar fraction depends on 
host galaxy properties such as Hubble type or B/T (Ode- 
wahn 1996; BJM08; Aguerri et al. 2009; M09; Weinzirl et 
al. 2009; Laurikainen et al. 2009), luminosity (Barazza et 
al. 2009; M09), stellar mass, and color (Nair & Abraham 
2010a; Cameron et al. 2010). Therefore we use comparison 
samples that are matched as well as possible to our Coma 
sample in Hubble type (SOs), luminosity (My < —18), 
color, and method of bar detection. We compare our re- 
sults for SOs to those of other studies in Coma and lower- 
density clusters (Abell 901/902 and Virgo). 

First, we compare to another study in the very dense 
environment of the central regions of the Coma cluster 
(galaxy number density n ~ 10,000 galaxies/Mpc 3 ) by 
T81. T81 uses visual classification on ground-based Kitt 
Peak National Observatory (KPNO) plates to detect bars 
in SOs brighter than My = —17.5 (very similar to our 
magnitude cutoff of My ~ —18). Therefore we compare 
his result to our optical bar fraction from visual classifica- 
tion (§ 13. 3p . Table [3] shows that our optical bar fraction 
from visual classification for SOs (60±11%) in the Coma 
core is higher than the result (42±7%) that T81 obtained 
after correcting raw galaxy counts for projection effect J^l. 
A clue to the reason for this difference comes from our find- 
ing that in many SOs in our sample, the bar ellipticity and 
its overall signature are diluted because the bulge is bright 
compared to the bar and it is large enough to encompass 
a large fraction of the bar length. In such cases, as dis- 
cussed in section § 13.2.21 the bar is harder to detect via 

and background objects contaminating the cluster field in the absence 
are spectroscopically-confirmed cluster members, we compare to the 
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any method, be it visual classification or ellipse fits, unless 
the image is of high quality and the classifier has signifi- 
cant expertise. In the case of T81, the visual classification 
was performed on ground-based optical plates, which are 
of lower quality than CCD images, making it even more 
difficult to detect such diluted or/and short bars. It is 
not possible to directly compare our case-by-case results 
with T81, as he does not publish the list of galaxies he 
classifies as barred and does not provide the lengths and 
ellipticities of the bars. However, we perform two indirect 
tests to gauge the impact of missing diluted, weak, and/or 
short bars. In our sample, seven of the 13 barred SOs have 
an observed peak bar ellipticity ebar < 0.4 (Fig. [T^]). If 
all of these galaxies were classified as unbarred, the bar 
fraction would drop to 30%. Alternatively, if the shorter 
(abar < 2 kpc) bars are excluded, the bar fraction would 
drop to 45%. These tests suggest that it is likely that the 
lower optical bar fraction of T81 is due to his missing some 
of these diluted or /and short bars. 

For a comparison to intermediate-density cluster envi- 
ronments, we use the study of M09 for the Abell 901/902 
cluster system (z ~ 0.165; n ~ 1000 galaxies/Mpc 3 , see 
Table 5 in Heiderman et al. 2009). To match our sample, 
we pick SO galaxies from the M09 study with My < —18, 
using classifications performed by the members of the 
STAGES collaboration for the Abell 901/902 cluster sys- 
tem (Gray et al. 2009; see Wolf et al. 2009 for more de- 
tails). In M09, inclined galaxies were picked as those with 
outer disk ellipticity edisk > 0.5, as traditionally done in 
bar studies using ellipse fitting. However, since we are 
only focusing on SO galaxies (which sometimes have large 
bulges/diffuse, thick stellar components that can dilute 
the edisk below 0.5 even for edge-on SOs, as discussed in 
§ I3.2.1[) we apply the same additional visual criteria out- 
lined in § l3.2.1l to the M09 SO sample to detect and remove 
highly inclined SOs. 

Figure [9] shows the host galaxy properties of the Abell 
901/902 and Coma SO samples. The two samples are well- 
matched in the mean luminosity, g — r color, and stel- 
lar mass. However, the results of a Kolmogorov-Smirnov 
(KS) test show differences in the overall distributions of 
g — r color and stellar mass (KS (p = 0.002, D — 0.4) and 
(0.005, 0.3), respectively), and are inconclusive for the dis- 
tributions of My (KS (p = 0.13, D = 0.2)). This is likely 
due to the fact that the Abell 901/902 sample has a tail of 
galaxies with masses both lower and higher than the Coma 
core sample, translating, respectively, into a tail of bluer 
colors and brighter absolute magnitudes. Oddly, there ap- 
pears to be a significantly larger offset between the Abell 
901/902 cluster and the Coma core if one uses the B — V 
color (Figure IHb) rather than the g — r color (Figure |Ht) : 
Abell 901/902 SOs appear ~ 0.2 mag bluer in B - V color 
compared to the Coma core, and the KS test suggests a 
large difference (KS (p = 6 * 10~ 16 , D = 0.8)). However, 
we believe that this B — V color offset is not real and is 
caused by the fact that the color transformations derived 
by Jester et al. (2005) for stars may not be adequate for 
SOs. 

In M09 bars were detected on optical {HST ACS 



F606W) images via the strict ellipse-fitting criteria only. 
We therefore also derive the bar fraction through visual 
classification and using the 'relaxed' ellipse-fitting criteria 
for the SOs in the M09 sub-sample in order to derive the 
corresponding bar fractions for comparison to Coma SOs. 
As shown in Tabled we find no statistically significant dif- 
ference in the SO bar fraction in Abell 901/902 and Coma 
clusters when detecting bars through visual classification, 
strict ellipse fit criteria, or relaxed ellipse fit criteria. For 
instance, from visual classification (Table [31 bottom sec- 
tion), the Coma bar fraction (60±11%) is slightly higher 
than the A901/902 bar fraction (55±5%), but the differ- 
ence is not statistically significant as the values are con- 
sistent within the error bars. Similarly, via the relaxed 
ellipse fitting criteria (Tabled middle section), the Coma 
bar fraction is 65±11%, while the A901/902 bar fraction is 
48±5%, (barely) consistent within the error bars. A com- 
parison of the bar and disk properties (such as ebar and 
Obar/^25) in Coma and the Abell 901/902 cluster system 
is discussed in § 13.71 

Next, we compare to results in Virgo from Erwin et al. 
(in preparation; Ell). Virgo is the most nearby cluster 
(D ~ 20 Mpc, z ~ 0.005) and is representative of a low- 
density cluster environment (n ~ 300 galaxies/Mpc 3 in 
the core region). We note however, that different environ- 
mental tracers paint different pictures in Virgo. While the 
number density (n ~ 300 galaxies/Mpc 3 ) is lower than 
that of Abell 901/902 (n ~ 1000 galaxies/Mpc 3 ) or the 
Coma core (n ~ 10,000 galaxies/Mpc 3 ), the velocity dis- 
persions in Virgo can be as high as 750 km/s (Binggeli 
et al. 1987), comparable to those seen in Abell 901/902 
and much higher than in groups (~ 100 km/s; Tago et 
al. 2008). These properties are relevant for the discus- 
sion of our results in § 13.81 Our Virgo comparison sam- 
ple consists of SO galaxies brighter than My = — 18 from 
Ell. Fig. [TOb shows that the two samples agree well in 
g - r color (KS (p = 0.3, D = 0.4)), but Fig. [lOp shows 
that the B — V colors of the Virgo SOs are significantly 
bluer (by ~ 0.15 mag, on average, and with KS values of 
(p = 8*10- n , D = 0.9)). We believe that this B-V color 
offset between Coma and Virgo is not real since the mea- 
sured g — r colors for the two samples agree well. We again 
think this offset is likely caused by the possibility that the 
color transformations derived by Jester et al. (2005) for 
stars may not be adequate for SOs. We compare optical 
bar fractions derived with all three methods: strict ellipse 
fitting criteria (ES), relaxed ellipse fitting criteria (ER), 
and visual classification performed by P.E., I.M., and S.J. 
according to the criteria outlined in § 13.31 Again, we do 
not find a statistically significant difference in the optical 
bar fraction (within the errors) for SOs in the Coma core 
and those in Virgo using any of the three bar-detection 
methods above (see Table l5]Fl 

A graphical representation of the trend of the bar frac- 
tion for SOs as a function of environment density is shown 
in Fig. QTJ We note that Fig. [TT] shows a hint of an in- 
crease in the mean bar fraction toward the dense core of 
the Coma cluster, however given the error bars, we cannot 
say whether this trend is significant. A comparison of the 



20 We note that Giordano et al. 2010 quote a much lower bar fraction (~ 30%) using visual classification for Virgo SO galaxies. This lower 
value is likely due to the fact that Giordano et al. 2010 include much fainter galaxies (down to Mb = —15), use a higher inclination cutoff 
(i = 73°), and a different method for selecting cluster members. 
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bar and disk properties, for SOs in Coma and the Virgo 
cluster is discussed in § 13.71 

We note that to compare to the lowest-density environ- 
ments (i.e., held galaxies), we would ideally like to have a 
comparison sample where bar detection is done quantita- 
tively via ellipse fits, and where the sample is matched 
to ours in both Hubble type (SO galaxies), luminosity 
(Mv ;$ —18), and color. However, there is as of yet no 
field comparison sample that fulfills all of the above re- 
quirements. The two large ellipse fit bar studies of field 
galaxies (BJM08 and Aguerri et al. 2009) are not adequate 
because they are mismatched in Hubble type and color 
(BJM08) or luminosity (Aguerri et al. 2009). Therefore, 
a comparison with these samples could be misleading, in 
light of recent results showing that the bar fraction varies 
non-monotonically with Hubble type, host galaxy lumi- 
nosity, and color (Nair & Abraham 2010a). A comparison 
of the bar fraction derived through visual classification to 
a matched subset of SOs from the RC3 is complex, because 
RC3 galaxies are a mix of field and Virgo cluster members. 
The best candidate for a field comparison is the recently 
released public catalog by Nair & Abraham 2010b, con- 
taining visual morphologies for ~ 14,000 SDSS galaxies. 
However, such a study is beyond the scope of this paper 
and therefore we defer this comparison to a later work. 

Recently Mendez-Abreu et al. (2010) used the Coma 
Treasury survey data to analyze the properties of barred 
galaxies in the Coma cluster (using visual classification to 
detect bars) . They do not select a disk sample but look for 
bars in all galaxies (including ellipticals and dwarfs). It is 
problematic and unconventional to quote the bar fraction 
from a sample of disk galaxies and ellipticals, particularly 
in the context of studying the bar fraction as a function 
of environment, as variations in this bar fraction can then 
be caused by the fact that the proportion of ellipticals to 
SOs to spirals changes strongly as a function of environ- 
ment. For this reason, our study and other studies quote 
the bar fraction as the fraction of disk galaxies hosting 
bars. Comparison of our work with the results of Mendez- 
Abreu et al. (2010) is therefore not straightforward, but 
nonetheless we attempt a comparison to check whether 
our findings are consistent. If we include all galaxies from 
our sample in the magnitude range —23 < M r < —14, to 
match their sample, we find a total (visual) bar fraction of 
7±2% (14/188), while Mendez-Abreu et al. (2010) find 9% 
(secure bars) and 14% (weak/uncertain bars). Although 
our results are broadly consistent within the uncertain- 
ties, there are several further caveats to this comparison. 
Mendez-Abreu et al. (2010) apply the axial ratio constraint 
b/a > 0.5 to their whole sample, regardless of morphology, 
while we only apply an inclination cutoff to our bright disk 
galaxies. In addition to the inclination cutoff, we also ex- 
clude SO galaxies deemed to be highly-inclined/edge-on by 
eye, since as discussed in § 13.2.11 using only the i > 60° 
cut misses highly-inclined and edge-on SOs with a more 
circular, thickened stellar component. Furthermore, in 
the process of selecting cluster members, Mendez-Abreu et 
al. (2010) apply a color cut where they discard all galaxies 
that have g — r color greater than 0.2 above their fit to the 
cluster red sequence. 

3.6. Observed and deprojected properties 



In Figure [T^] we show the observed (solid line) and de- 
projected (dotted line) distributions of the bar semi-major 
axis (abar) and bar ellipticity (ebar) for the 13 barred SO 
galaxies detected through ellipse fitting (using both strict 
and relaxed criteria). In MJ07, we found that deprojecting 
the bar semi-major axis and ellipticity for a large sample 
makes only a very small statistical difference (a factor of 
1.2), on average. However, since our Coma core sample 
is small, we deproject the observed radial profiles of ellip- 
ticity and PA and derive deprojected values of ebar and 
flbar- We perform the deprojection using a code developed 
by Laine et al. (2002) and used previously in Laine et 
al. (2002), Jogee et al. (2002a,b), and MJ07. 

We find an observed mean bar size of 2.5±1 kpc 
(2.9±1 kpc deprojected; Fig. [T3]) and the mean observed 
and deprojected bar ellipticity is 0.4±0.1. It is evident 
that deprojection docs not make a large difference in the 
mean abar and ebar- The observed and deprojected values 
of abar and ebar of the three bars detected through the 
relaxed criteria only are shown as filled and open circles, 
respectively in Fig. [T2] These galaxies satisfy the relaxed 
ellipse fit criteria, where the peak ellipticity of the bar 
is a local and not a global maximum due to the combi- 
nation of factors discussed in § 13.2.21 After deprojection 
removes the projection effects, which cause the ellipticity 
of the disk to be artificially boosted compared to the bar 
ellipticity, the peak bar ellipticity then becomes a global 
maximum in these three galaxies. Thus, after deprojec- 
tion, these three barred galaxies also pass the strict ellipse 
fit criteria. The observed and deprojected abar for these 
galaxies are shorter than average. We note, however, that 
the bars of the three galaxies that failed to meet the strict 
ellipse fit criteria before deprojection are not necessarily 
weak bars. The intrinsic ellipticities ebar of these three 
bars after deprojection are similar to the mean value of 
the whole sample (0.4±0.1). This can be understood from 
the combination of factors discussed in § 13.2.21 particularly 
the relative orientation of the bar with respect to the LON 
and projection effects. 

Fig. [T5] shows that the mean observed and deprojected 
bar lengths (abar) in our bright SO Coma core sample are 
similar to those of SOs in the Virgo cluster (Ell). The 
results from a KS test are consistent with similar distri- 
butions (giving KS {p = 0.5, D = 0.3) and (0.2, 0.4) for 
panels (a) and (b) of Fig[T31 respectively). We explore the 
comparison of the bar and disk properties as a function of 
environment density in more detail below. 

3.7. Properties of disks and bars in the Coma core 

We compare our sample of Coma cluster core SOs to 
the properties of SOs in less dense environments, namely 
to those in the intermediate-density Abell 901/902 clus- 
ter system (z ~ 0.165; M09) and those in the low density 
Virgo cluster (z - 0.005; Ell). 

Figure [T?] shows the distributions of (a) the galaxy 
disk i?25 (the isophotal radius where /ib reaches 
25 mag arcsec -2 ), (b) bar semi-major axis abar, mea- 
sured at the peak bar ellipticity ebar for all bars identi- 
fied through ellipse fitting (ER+ES), (c) the abar/^25 ra- 
tio, and (d) peak bar ellipticity ebar for the Coma, Abell 
901/902, and Virgo (Ell) SO samples. R 25 values for the 
Ell Virgo sample are from the RC3. R25 for the Coma 
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sample is estimated by ellipse fitting the galaxies on the 
ACS F475W images, which approximately correspond to 
SDSS g band. We calibrate the radial profiles of surface 
brightness to mag/arcsec 2 , then convert them from ~ g 
band (AB) to B mag arcsec -2 (Vega) using Equation 2. 
For two galaxies, we could not measure R25 radius be- 
cause a good fit could not be obtained of the outer disk of 
the galaxy due to the presence of a close companion. For 
the M09 Abell 901/902 sample, R 25 is calculated from the 
absolute Mb magnitudes according to 

log(-r^-) = -0.249 x M B - 4.00, (7) 
kpc 

from Schneider (2006). This formula is derived from an 
empirical relation measured for local spirals. To double- 
check its validity, we use it to calculate R25 radii for the 
Virgo SOs, where we already know R25 from RC3. Com- 
paring the calculated values with those from RC3 con- 
firms that the measured values from RC3 do follow the 
above relation, however it under-predicts the true R25 
by ~ 1.6 kpc on average. All three samples have sim- 
ilar mean bar and disk properties, but the bar semi- 
major axis and disk R25 distributions for Abell 901/902 
SOs have a tail to larger values. This tail corresponds 
to the tail of brighter SOs present in the Abell 901/902 
sample. The mean values of R25 for the Coma, Abell 
901/902, and Virgo SOs are 6.5±1.3 kpc, 5.9±2.8 kpc, 
and 5.0±2.0 kpc, respectively. The KS statistic reflects 
these differences in the Coma, STAGES, and Virgo distri- 
butions giving (p — 0.008, D = 0.5) between Coma and 
Virgo and (p = 10 , D — 0.4) between Coma and Abell 
901/902. The Abell 901/902 and Virgo SOs have similar 
mean abar/^25 ratios of ~ 0.4±0.16, although the range in 
values is large (~ 0.1-0.9). The KS statistic is consistent 
with similar distributions, giving (p = 0.9, D = 0.2) and 
(0.6, 0.3) between Coma and Abell 901/902 and Coma and 
Virgo, respectively. Coma SOs have a slightly lower mean 
abar/^25 = 0.35±0.12. An a har /R 25 ratio of - 0.3±0.2 
has also been found for field galaxies averaged over all Hub- 
ble types (e.g., MJ07, Menendez-Delmestre et al. 2007) 
and for SO galaxies (Erwin 2005). We note that the range 
of abar/ R25 spanned by the three samples is quite large 
(~ 0.1-0.9), however although our number statistics are 
small, this range is similar to that found for local field 
galaxies in MJ07. 

All three samples have very similar distributions in bar 
ellipticity (KS p ~ 0.9 and D ~ 0.2). The observed bar 
ellipticities we find for SOs in the Coma cluster as well as 
those for SOs in Virgo and Abell 901/902 are skewed to- 
ward lower values (e.g., mean ebar ~ 0.3-0.4) compared to 
the bar ellipticities in samples dominated by intermediate- 
to late-type galaxies (e.g., MJ07; BJM08; mean ebar ~ 0.5- 
0.7). This difference could be intrinsic (i.e., the bars in SO 
galaxies are really less elliptical than those in later Hubble 
types), or it could be due to the dilution by the bright 
bulges of the isophotes crossing the end of the bar, where 
the ellipticity is measured (see § I3.2.2[) . This effect has 
been demonstrated by Gadotti (2008). 

3.8. Discussion: implications for the evolution of SO bars 
and disks as a function of environment density 

What do our results imply for the evolution of bars and 
disks in bright SO galaxies as a function of environment? 



We first recapitulate our results. Using three detection 
methods (traditional ellipse fit criteria, relaxed ellipse fit 
criteria, and visual classification) , we found an optical bar 
fraction of 50±11%, 65±11%, and 60±11%, respectively 
for our sample of bright (My < —18) SO galaxies in the 
central region of the Coma cluster (§ 13. 4p . We find that the 
bar fraction and properties (e.g., ebar, «bar) in bright SO 
galaxies derived through all three of the above methods 
do not show a statistically significant variation (greater 
than a factor of ~ 1.3) between the dense central regions 
of Coma (n ~ 10,000 gal/Mpc 3 ), the intermediate-density 
Abell 901/902 clusters at z ~ 0.165 (n ~ 1000 gal/Mpc 3 ), 
and the low-density Virgo cluster (n ~ 300gal/Mpc 3 ; Ta- 
ble [3]) . We note that there is a hint that the mean bar 
fraction may show a slight increase as a function of envi- 
ronment density toward the dense core of the Coma clus- 
ter (Fig. Hip , however given the error bars, we cannot say 
whether this trend is significant. Below, we explore what 
our results may imply for the formation and evolution of 
bars. 

It has long been known that DM halo properties in- 
fluence bar formation and evolution. At high redshifts 
(e.g., z ~ 5-8), recent theoretical studies of galaxy evolu- 
tion using cosmological initial conditions find that bars are 
triggered by the triaxiality of DM halos and the asymmet- 
ric DM distribution as a whole (Romano-Diaz et al. 2008; 
Heller et al. 2007). These early bars are gas rich, and 
quickly decay. Subsequent bar generations form and are 
destroyed during the major- merger epoch (e.g., z ~ 2- 
4) due to the rapidly-changing potentials and gas dissi- 
pation associated with major mergers (Romano-Diaz et 
al. 2008). Although DM halos at early times can trigger 
bar formation due to their triaxiality, this triaxiality is di- 
luted as disks and other central components form. The 
DM halos become more symmetric, on a timescale that 
is a function of mass (e.g., Dubinski 1994; Kazantzidis et 
al. 2004; Heller et al. 2007). By z ~ 1 disks have also 
become more massive and stable. Simulations find that 
large-scale stellar bars forming at around this epoch are 
long-lived (Romano-Diaz et al. 2008; Heller et al. 2007). 
Interestingly, new observational results find that the bar 
fraction for the most massive disks (M* > IO^Mq) does 
not change between z ~ 0.6 and 0.2 (Cameron et al. 2010). 
However, the picture is complicated by the fact that for 
intermediate-mass disk galaxies (M* = 1O 1O 5 -1O 11 M0), 
the bar fraction builds up by a factor of two over that 
redshift range. In addition, at z ~ 0, the bar fraction and 
properties are a non-monotonic function of the host galaxy 
properties, such as stellar mass, luminosity, color, Hubble 
type, and SF history (BJM08; M09; Aguerri et al. 2009; 
Barazza et al. 2009; Weinzirl et al. 2009; Laurikainen et 
al. 2009; Gadotti 2011; Nair & Abraham 2010a). 

The picture above does not directly discuss environmen- 
tal effects. In fact, there are still few theoretical and ob- 
servational studies addressing this aspect of bar evolution. 
However, increasingly the emerging picture is suggesting 
that the frequency and properties of bars do not appear 
to be a sensitive function of environment (van den Bergh 
2002; Aguerri et al. 2009; M09; Barazza et al. 2009; 
Cameron et al. 2010; although see Giuricin et al. 1993 
and Elmegreen et al. 1990). 

How do the above results make sense in light of many 
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theoretical studies that show that galaxy interactions can 
trigger bars in unbarred galaxies (e.g., Noguchi 1988; Mi- 
hos & Hernquist 1996)? We present a tentative pic- 
ture below, considering the competing effects present in 
galaxy clusters. If a disk galaxy is sufficiently dynam- 
ically cold (i.e., Toomre Q < 1.5), it is susceptible to 
non-axisymmetric m = 2 instabilities (e.g., bars) whether 
spontaneously induced (e.g., Toomre 1981; Binney & 
Tremaine 1987) or tidally induced (e.g., Noguchi 1988; 
Hernquist 1989; Heller & Shlosman 1994; Mihos & Hern- 
quist 1996; Jogee 2006 and references therein). The ef- 
fect of the interaction depends on the geometry (i.e., pro- 
grade or retrograde encounter), with retrograde encoun- 
ters having little to no effect on an already existing bar 
(e.g., Gerin, Combes, & Athanassoula 1990; Steinmetz 
& Navarro 2002; Romano-Di'az et al. 2008; Aguerri & 
Gonzalez-Garci'a 2009). 

At z < 1.5, as clusters assemble and field galaxies fall 
into the existing cluster potential, let us now ask how the 
fraction and properties of bars in SOs might be expected to 
differ from the field environment. In a rich cluster, where 
the projected galaxy number density (n) and galaxy ve- 
locity dispersion (a) is high, the timescale for close inter- 
actions (or collision timescale, i C oii) will be short. We can 
estimate this timescale using: 



where n is the galaxy number density, <7 ga i is the galaxy 
velocity dispersion, and A is the cross-section for close in- 
teractions defined as 

A = 7r/(2r gal ) 2 . (9) 

We assume / is unity, r ga i ~ 10 kpc. For the Coma 
core CT ga i - 900 km/s, and n ~ 10,000 gal/Mpc 3 (The & 
White 1986) giving a short timescale for close interactions 
tooll ~ 90 Myr. 

However, although these close galaxy-galaxy interac- 
tions are frequent in a rich cluster, the large galaxy ve- 
locity dispersions present mean that each single encounter 
will be a high speed one. Unlike single slow, strong en- 
counters, a single high-speed encounter will typically not 
induce a large amount of tidal damage and not lead to 
major mergers. As a result, three factors may make it 
difficult for new bars to be induced in disk galaxies in a 
cluster. Firstly, single high-speed encounters may not be 
as effective in inducing bars as slow, strong encounters, 
because the timescale over which gravitational torques 
act is short. Secondly, over time, the cumulative effect 
of many high-speed and weak encounters (galaxy harass- 
ment), can tidally heat disks (e.g., Moore et al. 1996; 
Aguerri & Gonzalez-Garci'a 2009), making such disks dy- 
namically hot (with Toomre Q > 1.5), and thus less sus- 
ceptible to bar instabilities. Finally, in a cluster environ- 
ment, the accelerated star formation history (e.g., Balogh 
et al. 2004, Blanton et al. 2005b; Hogg et al. 2003) as 
well as physical processes such as ram pressure stripping 
(Gunn & Gott 1972; Larson et al. 1980; Quilis et al. 2000) 
will make SO disks gas-poor, thus making them less bar- 
unstable. We therefore speculate that these three factors, 
namely the predominance of high speed encounters over 
slow ones, the tidal heating of SO disks, and the low gas 
content of SOs in rich clusters, make it difficult for many 



new bars to be induced in SO disks as they infall from a 
field- like to a cluster- like environment. This scenario may 
explain, at least in part, our findings that there is no strong 
variation in the optical bar fraction of SOs across the range 
of low density to high density environments characterized 
by Coma, Virgo, and Abell 901/902 in our study, as well 
as claims by other studies that there is no difference in 
the bar fraction between clusters and the field (van den 
Bergh 2002, Aguerri et al. 2009, Barazza et al. 2009; M09). 
Our interpretation for the result that the bar fraction is 
not greatly enhanced in the dense Coma cluster is also in 
agreement with that of Mendez-Abreu et al. (2010). 

We note that it is possible that in rich clusters, the above 
effects, particularly the tidal heating, may cause existing 
bars to weaken. However, this effect is hard to robustly 
demonstrate observationally as the measured bar elliptic- 
ity is diluted by relatively large bulges in SOs (which dom- 
inate the disk population in clusters), while in the field, 
the disk population is dominated by spirals where such a 
dilution is not as severe (see ? !3.7p . 

It is also important to note that the arguments above, 
which explain why the bar fraction might not be greatly 
enhanced in rich clusters compared to the field, would lead 
to a rather different prediction for how the bar fraction in 
groups would compare to that in the field. In a group, the 
number density is moderately high (n ~ 10) but the galaxy 
velocity dispersions are typically low (a ~ 100; Tago et 
al. 2008). Therefore slow, strong encounters are expected 
to be frequent in groups. Such encounters are likely to in- 
duce extra bars in disk galaxies compared to the field, par- 
ticularly given the fact that the disks will not be stripped 
of their cold gas in groups as they would in rich clusters. 
In this context, we note that indeed higher bar fractions 
have been reported for early-type galaxies in binary pairs 
(Elmegreen et al. 1990) and early-type galaxies that are 
disturbed/interacting (Varela et al. 2004). 

4. BARS AND DISK FEATURES IN COMA DWARFS 

In addition to investigating bars in high-mass galaxies, 
we also take advantage of the exquisite resolution of the 
ACS (~ 50 pc at the distance of Coma) to search for bars 
and other disk features (e.g., spiral arms, edge-on disks) in 
the numerous dwarf galaxies in the central regions of the 
Coma cluster. Are some of these galaxies the remnants 
of late-type spirals that have gone through processing in 
a dense cluster environment? In Virgo, some early-type 
dwarfs are known to host features (e.g., lenses, bars, spiral 
arms; Sandage & Binggeli 1984; Binggeli & Cameron 1991; 
Jerjen et al. 2000; Barazza et al. 2002; Lisker et al. 2006; 
Lisker et al. 2007; Lisker & Fuchs 2009) suggesting the 
presence of a disk. Not only can such features provide 
clues to the formation history of these systems, but the 
presence or absence of bar structures has implications for 
the conditions necessary for bar formation and growth in 
galaxies (e.g., Mendez-Abreu et al. 2010). 

4.1. Identifying dwarf galaxies 

As outlined in §[51 we use a magnitude cut at Mi( 814 ) = 
— 18.5 (AB mag), (roughly corresponding to My = —18 
Vega mag) to separate dwarf and normal galaxies. A mon- 
tage of some of the faint, low-mass dwarfs in our Coma 
core sample is shown in Fig. 1151 Fig. 116a shows the dis- 
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tribution of absolute magnitude Mi( 814 ) of the galaxies in 
the faint dwarf sample. Fig. [TSb shows where the dwarf 
galaxies lie on a plot of the fi e (the surface brightness at 
R e ) vs. absolute magnitude Mj^gi^y Effective radii and 
fi c are from Hoyos et al. (2010), derived through single- 
component Sersic fits. 

Prior to applying the unsharp masking technique f§ !4.2[) 
on the dwarf sample, we pick out good candidates through 
a cut in surface brightness (/i c < 25 mag arcsec -2 ) and ra- 
dius (i?9o > 100 pc). These cuts remove very low surface- 
brightness objects, and those where we are unlikely to re- 
solve disk structure. We choose a size cut at Rgo > 100 pc, 
which is twice the ACS PSF at the distance of Coma. Out 
of the 417 dwarf galaxies, we find 333 dwarfs that satisfy 
these criteria. The fi c vs. Mj( 814 -) distribution of these 
galaxies is plotted with yellow triangles in Fig. IT5b . 

4.2. Identifying bars and other disk features in dwarfs 

In many dwarf galaxies disk features may not be readily 
apparent by eye (or traditional quantitative methods such 
as ellipse fitting) because their amplitude is very low and is 
overwhelmed by the smooth light from the galaxy. Which 
method is sensitive enough to detect faint spiral/bar struc- 
ture in such systems? Jerjen et al. (2000) used residuals 
from subtracting the azimuthally-averaged light profile of 
the galaxy from the original image to discover hidden spi- 
ral features in IC 3328 with deep VLT observations. They 
analyzed the spiral structure using Fourier expansion, find- 
ing that the amplitude of the spiral is only ~ 3 to 4%. 
However, upon further analysis of a larger sample of 20 
Virgo dwarfs, Barazza et al. (2002) find that some spiral 
or bar-resembling residuals may be artifacts from the com- 
bination of the increasing ellipticity and twisting isophotes 
(due to triaxiality) present in these galaxies and not ac- 
tual spiral structure. Fourier decomposition is similarly 
unsuccessful in many galaxies. Barazza et al. (2002) find 
that a much better method seems to be unsharp masking 
(e.g., Schweizer & Ford 1985; Mendez et al. 1989; Buta 
& Crocker 1993; Colbert et al. 2001; Erwin & Sparke 
2003). In this method, no assumptions about the light 
profile/inclination of the galaxy are necessary. Recently, 
Lisker et al. (2006) also successfully employed unsharp 
masking on ~ 470 Virgo dwarfs to look for evidence of 
bar/spiral structure. Graham et al. (2003) discovered two 
dwarf galaxies with spiral structure in the Coma cluster 
using unsharp masking as well as subtracting a symmetri- 
cal model to reveal non-symmetrical disk features (one of 
these galaxies is COMAil25937.988p28003.56 in Table IU 
while the other is not covered by the Coma ACS Treasury 
survey). Chilingarian et al. (2008) use unsharp masking to 
find disk features in dwarf galaxies in Abell 496. We there- 
fore use the unsharp masking method to seek out bar (or 
spiral) structures in the Coma cluster core dwarf sample. 

We perform unsharp masking for the 333 dwarf galaxies 
that fit the criteria outlined in ii l4.ll First, we smooth the 
galaxy images by convolving with a Gaussian using the 
IRAF task GAUSS. Then we divide the original galaxy 
image by the smoothed image. We choose the Gaussian 
smoothing kernel size to be ~ 25 pixels, corresponding to 
~ 625 pc for our galaxies. We also try a range of smooth- 
ing lengths from ~ 15-45 pixels (~ 375-1125 pc) for a 
subsample of the galaxies and find no substantial change 



in the results. A point made by Lisker et al. (2006) is that 
in some cases, it is desirable to use an elliptical smoothing 
aperture matched to the outer ellipticity and PA of the 
galaxy, in order to avoid spurious detections that resem- 
ble an edge-on disk. For this reason, in all cases where 
we suspect that the galaxy host an inclined/edge-on disk, 
we also perform the unsharp masking using an elliptical 
PSF to ensure that the structures found are not spurious 
detections. 

We find bars and/or spiral arms in 13 galaxies out of 
the 333 dwarfs in the unsharp-masked subsample. An ad- 
ditional eight galaxies show evidence of an inclined disk (or 
ambiguous inclined disk or bar). The galaxies where we 
find structure are listed in Table[4] Fig.[T7lshows examples 
of the residuals due to different types of structures: (a) spi- 
ral arms only, (b) bar and spiral arms, (c) inclined disk, 
(d) bar and/or spiral, (e) ambiguous bar or inclined disk 
embedded in a stellar halo. The galaxies with disk struc- 
ture are overplotted as cyan points in Fig.[16j Most (76%) 
of the galaxies where we find disk structure are brighter 
than M I ( 814 ) = — 16 (AB mag). We discuss the possible 
implications of these results below. 

4.3. Discussion: barred dwarf galaxies in the Coma core 

Using visual inspection and unsharp masking we find 
only 13 galaxies with bars and/or spiral arms in our 
Coma core dwarf subsample of 333 galaxies with /i c < 
25 mag arcsec -2 and Rqo > 100 pc. Does this result imply 
that faint/dwarf galaxies with disks are very rare within 
the Coma population, or rather that any existing disks in 
these galaxies are too dynamically hot to be unstable to 
disk instabilities? 

Studies have long been finding early-type dwarf galax- 
ies with spiral/bar structure in Virgo and Fornax (e.g., 
Sandage & Binggeli 1984; Binggeli & Cameron 1991; Jer- 
jen et al. 2000; Barazza et al. 2002; Lisker et al. 2006; 
Lisker et al. 2007; Lisker & Fuchs 2009). Lisker et 
al. (2006) search through 476 Virgo early-type dwarfs and 
find unambiguous stellar disk structure (bar/spiral) in 14 
of them, while another 27 have probable or possible disk 
features. Some authors have speculated that anywhere 
from 5% to 50% of Virgo early-type dwarfs have disk struc- 
ture, depending on the magnitude range under scrutiny 
(Lisker et al. 2006; Lisker & Fuchs 2009). 

Approaching the search for disks in early-type dwarfs a 
different way, Aguerri et al. (2005) investigate a sample 
of galaxies in Coma with —18 < Mb < —16 and classify 
them into two types dE or dSO depending on their sur- 
face brightness profile. Galaxies whose surface brightness 
profiles are well fitted by a single Sersic law are classified 
as dEs, and those with surface brightness profiles fitted 
with a Sersic plus exponential profile are classified as dSOs. 
Aguerri et al. (2005) find that about 30% of their Coma 
dwarf sample cannot be fitted well by a single Sersic law, 
suggesting that early-type dwarfs with disks may not be 
scarce in Coma. Graham & Guzman (2003) found evi- 
dence for outer disks in three out of a sample of 18 Coma 
early- type dwarfs, modeling the surface brightness profiles 
using a Sersic function in combination with either a cen- 
tral point source or a resolved central Gaussian component 
using high-resolution HST images. 

While it is still unclear whether all early-type dwarfs 
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with disk structure represent a distinct class of galaxies 
that are the product of a single formation mechanism, one 
plausible scenario is that they are formed through pro- 
cessing of faint, late-type spirals and irregulars in cluster 
environments (e.g., Kormendy 1985; Lin & Faber 1983; 
Graham et al. 2003; Lisker et al. 2006). This processing 
includes the loss of their gas through ram-pressure strip- 
ping (e.g., Gunn & Gott 1972) as well as the cumulative 
effects of harassment in a cluster environment (e.g., Moore 
et al. 1996). The simulations of Mastropietro et al. (2005) 
have shown that indeed late-type spirals can be repro- 
cessed into early-type dwarfs through cluster processes 
(such as harassment), and that these dwarfs do retain their 
stellar disk structure. On the other hand, observations 
of a small number of isolated early-type dwarf galaxies 
(e.g., Fuse et al. 2008; Hernandez-Toledo et al. 2010) ar- 
gue against cluster transformation processes as the sole 
explanation for the formation of these objects. 

We find evidence of only 13 dwarfs hosting disk instabil- 
ities (bar and/or spiral arms) in our unsharp-masked sub- 
sample of 333 dwarfs. This result is in broad agreement 
with the findings of Mendez-Abreu et al. (2010), who find 
a paucity of barred disks for Coma galaxies fainter than 
M r = —17. As also suggested by Mendez-Abreu et al. 
(2010), these results imply that although it is possible that 
as many as ~ 30% of dwarf galaxies in Coma may have a 
disk component (Aguerri et al. 2005), the majority do not 
have the necessary conditions to form or maintain bar and 
spiral instabilities, namely a disk that is dynamically cold. 
This result is consistent with previous studies, showing 
a paucity of thin disks in lower-luminosity dwarf galax- 
ies (Sanchez- Janssen et al. 2010; Yoachim & Dalcanton 
2006). 

5. SUMMARY 

We use ACS F814W images from the Hubble Space Tele- 
scope ACS Treasury survey of the Coma cluster at z ~ 0.02 
to study the fraction and properties of barred galaxies in 
the central region of Coma, the densest environment in the 
nearby Universe. The available data span 274 arcmin 2 , 
where approximately 75% of the data are within 0.5 Mpc 
of the cluster center, and contain thousands of sources 
down to a limiting magnitude of I = 26.8 mag in F814W 
(AB mag). We initially select 469 cluster members and 
split the sample with a magnitude cut at My < —18 
(Vega mag). Using this magnitude cut, we investigate 
two different regimes: (1) the fraction and properties of 
bright My < —18 SO galaxies and (2) the presence of 
bars and other disk features (e.g., bars and spiral arms) in 
faint / dwarf galaxies in the Coma core. Our results for the 
two populations are described below. 

(1) For SO galaxies: We select a sample of 32 bright 
SO galaxies based on visual classification supplemented 
by multi-component decompositions in ambiguous cases 
(§ 13. ip . After discarding 12 highly-inclined galaxies, 
we identify and characterize bars in the remaining 20 
moderately-inclined SOs using three methods: ellipse fits 
where the bar is detected through strict criteria (the peak 
bar ellipticity ebar is required to be a global maximum in 
the radial profile of ellipticity); ellipse fits where the bar 
is detected through relaxed criteria (which do not require 
the peak bar ellipticity ebar to be a global maximum); and 



visual classification. We find: 

• The optical bar fraction for our bright SO sample is: 
50± 11%, 65± 11%, and 60± 11% based on ellipse 
fits with traditional and relaxed criteria, and visual 
classification, respectively (Table [5]). 

• We compare to results from studies in less dense 
environments (Abell 901/902 and Virgo) and find 
that the bar fraction, as well as the mean quan- 
titative properties of the SO bars and disks (e.g., 
R25, flbar, ebar) do not show a statistically signifi- 
cant variation, within the error bars, for samples of 
matched SOs in environment densities ranging from 
n ~ 300 gal/Mpc 3 (Virgo), n ~ 1000 gal/Mpc 3 
(Abell 901/902), and n ~ 10,000 gal/Mpc 3 (Coma), 
with high galaxy velocity dispersions a ~ 800 km/s 
(Table H Fig. EJ}. We note that there is a hint 
that the mean bar fraction may show a slight in- 
crease as a function of environment density toward 
the dense core of the Coma cluster (Fig. [Til , how- 
ever given the error bars, we cannot say whether 
this trend is significant. We speculate that the bar 
fraction among SOs is not dramatically enhanced in 
rich clusters compared to low density environments 
due to several factors. Firstly, SOs in rich clusters 
are likely to be more stable to bar instabilities be- 
cause they are dynamically heated by the cumu- 
lative effect of many high-speed, weak encounters 
(galaxy harassment), and additionally are gas poor 
as a result of ram pressure stripping and accelerated 
star formation. Secondly, individual high-speed en- 
counters in rich clusters may be less effective than 
individual slow strong encounters in inducing bars. 
The combination of these effects precludes an en- 
hancement in the bar fraction for SO galaxies in 
cluster environments compared to the field. Our 
results are in agreement with recent observational 
studies which find no difference in the fraction of 
barred galaxies with environment density over all 
Hubble types. 

(2) For faint/dwarf galaxies: We select a sample of 
417 galaxies fainter than Mt( 814 ) = —18.5 (AB mag; 
§ 14.11) where we utilize our ~ 50 pc resolution to look 
for disk structures such as bars and spiral arms using 
visual classification of unsharp-masked images. After ap- 
plying unsharp masking to a subsample of 333 dwarfs 
(/x e < 25 mag arcsec -2 , Rgo > 100 pc; § I4.2[) . we find only 
13 dwarf galaxies with a bar and/or spiral arms, and an 
additional eight galaxies where an inclined disk may be 
present (Fig. ITT]) . These results suggest that either disks 
are not common in these galaxies in the Coma cluster core, 
or that any disks present are too hot to form instabilities. 
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Table 1 

Optically barred bright (M v < —18) SOs. 



Bar ID 


Galaxy ID 


RA 


DEC 


^1(814) Bar detection method 


ebar 


ebar 


flbar 


Obar 






(AB mag) 




obs 


dcp 


obs (kpc) 


dep (kpc) 


(1) 


(2) 


(3) 


(4) 


(5) (6) 


(7) 


(8) 


(9) 


(10) 



Bl 


COMAil25710.767p272417.44 


194.29486 


27.404846 


-20.1 


ES, ER, V 


0.35 


0.26 


2.97 


2.99 


B2 


COMAil25833.136p272151.77 


194.63806 


27.364383 


-20.0 


ES, ER, V 


0.68 


0.54 


3.19 


3.52 


B3 


COMAil25928.728p28225.90 


194.86970 


28.040531 


-19.7 


ES, ER 


0.49 


0.40 


2.26 


2.44 


B4 


COMAil25929.956p275723.26 


194.87481 


27.956462 


-21.1 


ES, ER, V 


0.47 


0.51 


4.36 


5.04 


B5 


COMAil25946.775p275825.88 


194.94490 


27.973856 


-20.6 


ES, ER, V 


0.39 


0.70 


2.22 


2.82 


B6 


COMAil25956.707p275548.62 


194.98628 


27.930173 


-19.9 


ES, ER, V 


0.55 


0.46 


3.42 


4.00 


B7 


COMAil3022.156p28249.08 


195.09231 


28.046968 


-20.6 


ES, ER, V 


0.37 


0.46 


1.92 


2.23 


B8 


COMAil3038.731p28052.22 


195.16137 


28.014507 


-20.3 


ES, ER, V 


0.47 


0.31 


3.57 


3.77 


B9 


COMAil3042.753p275816.88 


195.17814 


27.971355 


-21.3 


ES, ER, V 


0.25 


0.31 


1.46 


1.55 


B10 


COMAil3042.833p275746.98 


195.17846 


27.963052 


-20.3 


ES, ER, V 


0.43 


0.46 


2.48 


2.64 


Bll 


COMAil25930.825p275303.42 


194.87843 


27.884283 


-21.1 


ER, V 


0.33 


0.35 


1.99 


2.33 


B12 


COMAil3017.020p28350.10 


195.07092 


28.063917 


-19.0 


ER, V 


0.29 


0.40 


1.14 


1.45 


B13 


COMAil3027.971p275721.54 


195.11654 


27.955985 


-20.1 


ER, V 


0.31 


0.47 


2.09 


2.79 



Note. — (1) Bar ID; (2) Galaxy ID as given in the Coma Treasury survey DR2 (Hammer et al. 2010); (3) RA (J2000); (4) DEC (J2000); (5) M I(8 i 4) 
absolute magnitude in AB mag; (6) Bar det ection method: 'ES' - strict ellipse fit criteria (§ 13.2.21 , 'ER' - relaxed ellipse fit criteria f§ I3.2.2K and 'V 
- visual classification on direct image f ^ 13. 31 (7) Observed peak bar cl lipti city ebar; (8) Dcprojcctcd peak bar cllipticity; (9) Observed bar semi-major 
axis a-bar measured at ebar! (10) Dcprojcctcd bar semi-major axis ("§ 13.711 ■ 
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Table 2 

Optical bar fraction for bright (My < —18) SOs based on different methods. 



Method 


Unbarred 


Barred 


/bar, opt 


Ellipse fit, strict 


10 


10 


50±11% 


Ellipse fit, relaxed 


7 


13 


65±11% 


Visual classification 


8 


12 


60±11% 



Note. — optical bar fraction for the 20 moderately in- 
clined (i < 60°) bright (My < -18) SO galaxies. Barred 
galaxies arc characterized through: (1) ellipse fitting using 
the strict criteria (where e^ar ls required to be a global 
maximum in the ellipticity profile), (2) ellipse fitting us- 
ing relaxed criteria (ebar can be a local maximum), and 
(3) visual classification. 



Table 3 

Optical bar fraction for bright (My < —18) SOs in different environments. 

Study Environment Number density (gal/Mpc 3 ) Velocity dispersion (km/s) SO /bar for My < —18 

(1) (2) (3) (4) (5) 



Bars identified through ellipse fitting (strict criteria; § I3.2.2[) 



this work 


Coma central region, z ~ 0.02 


10,000 


900 (o) 


50±11% (10/20) 


M09 (b) 


Abell 901/902 clusters, z ~ 0.165 


1000 


800-1200 (c> 


39±5% (38/98) 


Ell 


Virgo, z ~ 


300 


400-750 (<1> 


44±14% (12/27) 




Bars identified through ellipse fitting 


; (relaxed criteria; 


S 13.2.21 




this work 


Coma central region z ~ 0.02 


10,000 


900 


65±11% (13/20) 


M09 (i>) 


Abell 901/902 clusters, z ~ 0.165 


1000 


800-1200 


48±5% (47/98) 


Ell 


Virgo, z ~ 


300 


400-750 


48±14% (13/27) 



Bars identified through visual classification 



this work Coma central region, z ~ 0.02 10,000 900 60±11% (12/20) 

T81 (e) Coma central region, z ~ 0.02 10,000 900 42±7% (19/45) 

M09 (6) Abell 901/902 clusters, z ~ 0.165 1000 800-1200 55±5% (54/98) 

Ell (/) Virgo, 2 ~ 300 400-750 59±9% (16/27) 



Note. — T81: Thompson 1981; M09: Marinova et al. 2009; Ell: Erwin et al. (in prep.) 

(a) The & White (1986) 

(b) We use a sub-sample from M09, with the criteria outlined in § 13.51 

(c) Heiderman et al. (2009) 

(d) Binggeli et al. (1987) 

(c) Bar classification is performed on ground-based KPNO plates. 

(f) For this paper, visual classification is performed on the Ell sample by P.E., I.M., and S.J. using the criteria outlined in § 13.31 
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Table 4 

Galaxies in the faint (My > —18) sample where we find disk structure through unsharp masking. 



Galaxy ID 


RA 


DEC 


Visit 


•^1(814) 


Structure 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 



COMAil3030.949p28630.18 

COMAil3041.192p28242.38 

COMAil3018.883p28033.55 

COMAil3013.398p28311.81 

COMAil3007.727p28051.91 

COMAil25930.062p28237.71 

COMAil25904.792p28301.21 

COMAil25911.545p28033.38 

COMAil25953.930p275813.76 

COMAil25937.988p28003.56 

COMAil3035.418p275634.05 

COMAil3024.823p275535.89 

COMAil25950.181p275445.54 

COMAil25820.533p272546.03 

COMAil25815.275p272752.96 

COMAil25814.969p272744.81 

COMAil25825.308p271200.04 

COMAil25623.788p271402.30 

COMAil25638.099p271304.09 

COMAil25845.297p274650.75 

COMAil25845.906p274655.90 



195 


12895 


28. 


10838 


02 


195. 


17163 


28. 


04510 


08 


195. 


07867 


28. 


00931 


09 


195. 


05583 


28. 


05328 


10 


195. 


03219 


28. 


01441 


10 


194. 


87525 


28. 


04380 


13 


194 


76997 


28. 


05033 


14 


194 


79811 


28. 


00927 


14 


194 


.97471 


27. 


97048 


18 


194. 


90827 


28. 


00098 


19 


195. 


14758 


27. 


94279 


22 


195. 


10342 


27. 


92663 


23 


194 


.95909 


27. 


91265 


25 


194 


.58555 


27. 


42945 


45 


194 


.56364 


27. 


46471 


45 


194 


.56237 


27. 


46244 


45 


194 


.60545 


27. 


20001 


59 


194 


.09912 


27. 


23397 


63 


194 


.15875 


27. 


21780 


63 


194 


.68873 


27. 


78076 


75 


194 


.69126 


27. 


78219 


75 



-17.6 


inclined disk 


-17.8 


spiral 


-18.1 


bar 


-17.7 


inclined disk 


-13.3 


inclined disk 


-13.4 


bar/inclined disk 


-18.3 


bar 


-18.2 


inclined disk 


-16.9 


spiral 


-18.0 


bar -f- spiral 


-17.7 


bar 


-18.2 


bar -f- spiral 


-17.5 


bar 


-18.3 


spiral 


-17.2 


inclined disk 


-15.5 


bar -f- spiral 


-18.4 


bar 4- spiral 


-18.1 


bar 4- spiral 


-16.4 


bar 


-15.3 


bar/inclined disk 


-14.7 


bar/inclined disk 



Note. — (1) Galaxy ID as given in the Coma Treasury survey DR2 (Hammer et al. 2010); 
(2) RA (J2000); (3) DEC (J2000); (4) HST visit number; (5) M I(8 i 4) absolute magnitude in AB 
mag; (6) type of disk structure detected through unsharp masking (see § |4.2[> . 
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Fig. 1. — (a) Absolute magnitude (Mjfgj^) distribution of the bright, non-dwarf (black solid line; My < —18) and dwarf (green dash-dot 
line; My > —18) galaxies in our Coma core cluster member sample (§EJ- Most galaxies are dwarfs with My > —18. (b) Central surface 
brightness /io (from the S Extractor source catalogs in DR2) vs. absolute magnitude Mjfgi4) for the bright, non-dwarf (black circles) and 
dwarf (green plus) cluster core samples. 
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Fig. 2. — For the 10 visually-ambiguous E/SO galaxies, single and multiple component decompositions were performed with GALFIT on 
the 2D light distribution to help determine whether they are an elliptical (E) or SO. This figure shows examples of the data (left), GALFIT 
model (middle), and (data-model) residuals (right) for three of these galaxies, which ended up being classified after decomposition (§CTT} as 
E (top row), SO/a (middle row), and still ambiguous E/SO (bottom row). For the galaxy (COMAil25930.268p28115.17) in the top row the 
best fit model consists of a Sersic component of half light radius R e ~ 1.4 kpc, within which is embedded a compact disk and a nuclear point 
source. This galaxy does not make the cut to be an SO as it lacks a bulge and extended outer disk, and it is not included in the analysis 
of the bar fraction for SOs. The galaxy in the second row (COMAil25938.323p275913.84) is classified as an SO/a because after fitting a 
(Sersic+Exponential) model, corresponding to a bulge and extended outer disk, one can see extended spiral structure in the residuals. Since 
spiral structures only exist in disks, this confirms the presence of an outer disk. We classify this galaxy as an SO/a rather than an Sa because 
the spiral structure is revealed in the residuals and is not readily visible on the direct image. The galaxy (COMAil25950.103p275529.47) 
in the third row is classified as an ambiguous E/SO because, based on the residuals and other factors (x 2 and fit parameters), it is still not 
possible to determine whether this galaxy is more likely to be an elliptical with an inner debris disk or an SO. 
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Fig. 3. — (a) Absolute magnitude distribution of bright (Mv < —18) galaxies in our Coma core sample (§ 13. It . (b) Distribution 
of morphological types (E, E/SO, SO, Sp) in the bright non-dwarf sample. SO galaxies (32), comprise 94% of the bright disk galaxies. 
Morphological types are from visual classification, supplemented with 2D, multi-component decomposition for visually ambiguous cases 
(ij 13.111 . (c) Stellar mass distribution of the bright galaxies in our Coma core sample, with the SO galaxies shown in green. The SOs have 
masses between 10 9 ' 5 and 10 11 Mq. (d) g — r color-magnitude diagram of the bright cluster sample and a subset (30%) of the dwarf sample 
with available SDSS magnitudes. We overplot the relation from Blanton et al. (2005) for the break between the red sequence and blue cloud. 
Most elliptical and SO galaxies lie on the red sequence. 
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Fig. 4. — Examples of some SO galaxies from our bright sample (§ 13. 11 . The scale bars show 1 kpc. The galaxies shown are: (a) CO- 
MAil25704.336p273133.26, (b) COMAil25710.767p272417.44, (c) COMAil25833.136p272151.77, (d) COMAil25832.060p272722.85, (e) CO- 
MAil25928.728p28225.90, (f) COMAil25929.404p275100.51, (g) COMAil25929.956p275723.26, (h) COMAil25930.825p275303.42, (i) CO- 
MAil25931.455p28247.62, (j) COMAil25932.789p275900.95, (k) COMAil25938.323p275913.84, (1) COMAil25939.657p275713.86. 
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Stfl - 




Fig. 5. — Left: HST image and radial profiles of surface brightness, e, and PA of a highly-inclined SO with an outer, diffuse, thickened 
stellar component (see § 13. 2. Il l, Right: the ellipse fits are overlaid onto the galaxy image. The top two panels are shown with a stretch that 
enhances the thin disk and boxy bulge, while the bottom panel shows the outer disk. The thickened, diffuse, outer stellar component causes 
the outermost isophotes to have e ~ 0.4, which is less than the quantitative in clinat ion cut of e > 0.5. Therefore, we classify this galaxy as 
highly-inclined using visual classification according to the criteria outlined in § 13.2.11 
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Fig. 6. — Left: HST image and radial profiles of surface brightness, e, and PA of a barred cluster galaxy. In this example, the traditional 
bar signature is evident in the smooth rise of the e to a global maximum of ~ 0.4, while the PA remains relatively constant in the bar region. 
The e then drops and the PA changes, indicating the transition to the disk region. Right: the ellipse fits are overlaid onto the galaxy image. 
The top two panels are shown with a stretch that enhances the inner disk and bar regions, while the bottom panel shows the outer disk. See 
§1323] for details. 
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Fig. 7. — The panels are as in Fig. [6] but here we show an example of a barred galaxy that does not meet the strict ellipse-fit criterion 
requiring that e bar is the global maximum in the e radial profile. In this case, the observed outer disk cllipticity e^k is higher than e bar , 
making it a local maximum in the e radial profile. This happens due to a combination of properties of the galaxy: (1) the galaxy is inclined 
(i ~ 51°) causing the outer disk to be elongated along the line of nodes with a significant ellipticity (ejisk = 0.37); (2) the stellar bar is 
significantly offset (by ~ 45°) with respect to the line of nodes and hence its intrinsic axial ratio is diluted by projection effects; (3) the stellar 
bar has a significant fraction of its length inside a very luminous bulge, and the measured bar ellipticity is diluted to lower values than the 
true ebar- Therefore this galaxy is identified as 'barred' through the relaxed ellipse- fitting criteria. We find three such cases among the bright 
Coma SO galaxies (see § 13. 2. 211 . 



24 




Fig. 8. — Barred bright {My — 18) SO galaxies in the Coma cluster core found through ellipse fits (strict and relaxed criteria) and visual 
classification (see £j I3,2l and Table [TJ. All bars identified through the relaxed ellipse fit criteria are also identified by visual classification and 
vice versa. Bright stars such as the one in B8 are masked during the fitting. 
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Fig. 9. — The properties of bright SO galaxies in the Abell 901/902 cluster system (dashed pink lines) and Coma (solid black line). The 
vertical lines show the mean values for each distribution. The two samples are well-matched in mean luminosity, g — r color, and stellar 
mass. However, the results of a Kolmogorov-Smirnov (KS) test show differences in the overall distributions for g — r color and stellar mass 
(KS (p = 0.002, D = 0.4) and (0.005,0.3), respectively) and are inconclusive for the distributions of My (KS (p = 0.13, D = 0.2)). This is 
likely due to the fact that the Abell 901/902 sample has a tail of galaxies with masses both lower and higher than the Coma core sample, 
translating, respectively, into a tail of bluer colors and brighter absolute magnitudes. The Abell 901/902 SOs appear ~ 0.2 mag bluer in B — V 
color, on average (KS (p = 6 * 10 -16 , D = 0.8)). We believe that this B — V color offset is not real and is caused by the fact that the color 
transformations derived by Jester et al. (2005) for stars may not be adequate for SOs. 
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Fig. 10. — The properties of bright SO galaxies in the Virgo cluster (dashed green line) and Coma (solid black line). The vertical lines show 
the mean values for each distribution. The two samples are well-matched in mean luminosity and g — r color. A KS test shows differences 
in the overall distribution of Mb absolute magnitude (p = 0.02, D = 0.4), whereas the distributions of g — r color are similar (KS p = 0.3, 
D = 0.3). While the Virgo SOs appear bluer in B — V (~ 0.15 mag on average, with KS (p = 8 * 10~ n , D = 0.9)) we believe that this color 
offset is not real since the measured g — r colors (from SDSS in panel (c)) for the Coma and Virgo samples agree well. A possible reason for 
the offset in B — V is that the the color transformations derived by Jester et al. (2005) for stars may not be adequate for SOs. 
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Fig. 11. — The optical bar fraction for SO galaxies characterized through three methods (ellipse fitting with strict criteria, ellipse fitting 
with relaxed criteria, and visual classification) as a function of environment density. The different environments probed are the high-density 
core of Coma (n ~ 10,000 gal/Mpc 3 ), the intermediate-density Abell 901/902 cluster system (n ~ 1000 gal/Mpc 3 ), and the low-density Virgo 
cluster (n ~ 300 gal/Mpc 3 ; § 13.51 1. The bar fraction for SOs does not show a statistically significant variation across the environments probed, 
within the error bars. 
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Fig. 12. — Observed (solid line) and deprojected (dotted line) bar size (a^^) (a) and ellipticity e^ar (b) distributions for the 13 barred 
SO galaxies detected through ellipse fitting (including ones detected through relaxed criteria). The observed and deprojected values for the 
three bars detected through the ellipse fit relaxed criteria are shown as filled and open circles, respectively (§ I3,2.2j l, The mean observed 
a bar for our barred SOs (including those detected with relaxed criteria) is 2.5±1 kpc (2.9±1 kpc deprojected), while the mean observed and 
deprojected et, ar is 0.4± 0.1. Most (85%) of bars have an observed et, ar < 0.5. We also note that all extra bars that were detected via the 
relaxed ellipse fit criteria on the observed images, would be detected via the strict ellipse fit criteria after deprojection. This is due to the 
fact that the latter removes projection effects, which cause the maximum bar ellipticity eb ar to go from a local maximum in the radial profile 
of ellipticity to a global maximum. 
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Fig. 13. — Comparison of the observed (a) and deprojected (b) bar semi-major axis ab ar distributions for our sample of barred bright 
Coma core SOs and those in Virgo from Erwin et al. (in prep.). For both samples, the distributions include barred galaxies detected through 
the strict ellipse-fitting criteria as well as the relaxed criteria (§ I3.2.2H . The vertical lines show the mean values for each distribution. We do 
not find a significant difference in the mean observed and deprojected bar size between Coma and Virgo SOs. The results from a KS test are 
consistent with similar distributions (giving KS (p = 0.5, D = 0.3) and (0.2, 0.4) for panels (a) and (b), respectively. 
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Fig. 14. — Distributions of (a) disk R25 (the isophotal radius where /xb reaches 25 mag arcsec - 2 ), (b) bar semi-major axis measured 
at the peak bar ellipticity et, ar for all bars identified through ellipse fitting (ER+ES), (c) at, ar /-R25 ratio, and (d) peak bar ellipticity eb ar 
for the Coma SO sample (solid black) and the comparison samples of SOs from the intermediate-density cluster system Abell 901/902 (dashed 
pink) and the low-density Virgo cluster (dotted green). R25 values for Coma and Abell 901/902 SOs are derived as described in § 13.71 while 
R25 for Virgo galaxies are from the RC3. The vertical lines show the mean values for each distribution. All three samples have similar mean 
bar and disk properties, but the bar semi-major axis and disk R25 distributions for Abell 901/902 SOs have a tail to larger values. The KS 
statistic reflects the differences in the Coma, STAGES, and Virgo distributions of R25, giving (p = 0.008, D = 0.5) between Coma and Virgo 
and (p = 10" 4 , D = 0.4) between Coma and Abell 901/902. 




Fig. 15. — Examples of some of the faint, low-mass dwarfs in the Coma core sample. The scale bars are 1 kpc. The galaxies 
shown arc: (a) COMAil3011.143p28354.92, (b) COMAil3025.977p28344.68, (c) COMAil3026.152p28032.02, (d) COMAil3029.853p28400.85, 
(c) COMAil3030.027p28135.08, (f) COMAil3039.068p28437.52, (g) COMAil3041.192p28242.38, (h) COMAil3047.670p28533.95, (i) CO- 
MAil3048.045p28557.42, (j) COMAil3050.590p28356.56, (k) COMAil3052.942p28435.86, (1) COMAil3030.949p28630.18. 
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Fig. 16. — (a) Absolute magnitude (iWi(gi4)) distribution and (b) plot of surface brightness vs. absolute magnitude (Mi( gl 4)) of the 417 
galaxies in the Coma faint sample. The cyan points show the values for the 21 faint galaxies where we find disk structure (bar, spiral, edge-on 
disk) through unsharp masking. Most (76%) of the objects where we find disk structure have Afj/gj^j < —16. 
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Fig. 17. — Examples of five galaxies from the dwarf (My > —18) sample highlighting the different types of disk structure that we find 
through unsharp masking: (a) spiral arms, (b) bar+spiral arms, (c) edge-on disks, (d) bar and/or spiral structure, and (e) ambiguous 
bar/edge-on disk (§ 14.20 . We use a Gaussian smoothing kernel size of ~ 25 pixels, corresponding to ~ 625 pc at the distance of Coma. The 
original HST images are shown in the left panels and the corresponding residuals highlighting the disk structure are on the right. The scale 
bars show 1 kpc. 



